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Abstract 

In the node-weighted prize- collecting Steiner tree problem (NW-PCST) we are given an undi- 
rected graph G = {V, E), non-negative costs c{v) and penalties tt{v) for each v £V. The goal is 
to find a tree T that minimizes the total cost of the vertices spanned by T plus the total penalty 
of vertices not in T. This problem is well-known to be set-cover hard to approximate. Moss 
and Rabani (STOC'Ol) presented a primal-dual Lagrangean-multiplier-preserving 0(ln |F|)-ap- 
proximation algorithm for this problem. We show a serious problem with the algorithm, and 
present a new, fundamentally different primal-dual method achieving the same performance 
guarantee. Our algorithm introduces several novel features to the primal-dual method that may 
be of independent interest. 

1 Introduction 

The node-weighted Steiner tree problem (NWST) is a fundamental and well understood network 
design problem, where we are given an n-node, undirected graph G = iy,E), a non- negative cost 
c{v) for each vertex v ^ and a set of terminals R QV . The goal is then to find a tree T D R 
that has minimum total cost X]^gj'C(f). 

In this paper, we consider the prize- collecting version of the problem. As in standard NWST, 
the input in an instance of the node-weighted prize- collecting Steiner tree problem (NW-PCST) 
again consists of an n-node, undirected graph G = (V, E) and a non-negative cost c{v) for each 
vertex v , but instead of a set of terminals we have now a non-negative penalty value 7r(f ) for 
each node v . We wish to find a tree T that minimizes 

^c{v)+ ^ tt{v). 

veT v&V\T 

Both NWST and NW-PCST have numerous practical and theoretical applications (e.g., j5l|9l 
116]). and are well-known to be NP-hard. From an approximation point of view, there is a relatively 
straight-forward, approximation-factor preserving reduction from the sei-couer problem, and there- 
fore no o(ln n)-approximation algorithm exists for either one, unless NP C DTIME(nP°^y^°S(")) 
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[6l [TT]. In fact, NWST and NW-PCST are significantly harder than the corresponding edge- 
weighted variants, for which constant approximation algorithms are known jJl \TU[ \TE[ El [3l [1]. 

Klein and Ravi [TT| showed that NWST admits an O(lnn) approximation, matching the above 
in-approximability bound. Guha et al. [9] later gave a primal-dual interpretation of this algorithm 
using a natural linear programming formulation of the problem. 

The focus of this paper will be the approximability of the prize-collecting variant of the problem. 
We will in fact address the rooted version of the problem, where a specific root vertex r £ V has to 
be part of the output tree T. Clearly, any algorithm for the rooted NW-PCST immediately yields 
an algorithm with the same performance guarantee for the unrooted version. The main result of 
this paper is stated in the following theorem. 

Theorem 1.1. There is a Lagrangean multiplier preserving (LMP) O (hi n)- approximation algo- 
rithm for NW-PCST; i.e., there is an algorithm that, given an NW-PCST instance, computes a 
tree T containing the root r such that 



where opt is the value of an optimum solution to the problem, and a = O(lnn). 

We note that it is reasonably straight forward to obtain a non-LMP 0(lnn)-approximation al- 
gorithm for NW-PCST via a standard threshold rounding approach for the natural LP formulation 
of the problem [Tj. Ensuring the LMP property is highly non-trivial, however, and of crucial impor- 
tance in the design of approximation algorithms for partial versions of NWST via the Lagrangean 
framework of Chudak et al. [4j. Two such partial problems that were considered in [T3] are the 
quota and budget versions of NWST. In the former, we are given a non-negative profit p{v) for each 
vertex v € V and a quota Q > 0, and we wish to find a connected set T of smallest cost whose 
vertices have profit at least Q. In the latter problem, we are given a budget B > 0, and wish to 
find a connected set T of largest total profit whose cost is at most B. 

The correctness of Theorem |1.1| was previously claimed by Moss and Rabani |13l I14j : the 
primal-dual algorithm presented there does, however, have a crucial technical mistake that appears 
not to have a simple fix, as we explain in the next subsection. For this reason, we present a 
fundamentally different primal-dual approach. Beside proving correctness of Theorem and 
therefore establishing correctness of results that rely on using Theorem |1.1| (e.g. approximation 
results for the previous mentioned partial NWST problems [H] or for problems arising in the study 
of contagion processes in networks |12tl8j), our algorithm features several new ideas that it adds to 
the known primal-dual repertoire, and therefore might be of independent interest. 

Wc begin by presenting the natural integer linear programming formulation for NW-PCST (see 
also p^). Afterwards, we give a high-level description of Moss and Rabani's algorithm, and show 
an example where it does not perform correctly. We then outline our algorithm, highlighting the 
main novelties with respect to currently known primal-dual approaches. 

1.1 LP Formulation for NW-PCST 

The integer program has a variable for each vertex v €z V that has value 1 if u is part of the 
output tree T, and x^ = otherwise. We also have variables for all S" C y where V = V \ {r}. 
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We let ^5 = 1 if S" is the set of vertices not spanned by T, and zs = otherwise. In the fohowing 
we let T{S) for a set 5 C 1/ be the set of all vertices u £ V \ S that have a neighbour in S, and we 
let vr(5) be equal to J2ves 



min ^ c{v)xy + ^ TTiS)zs (P) 
vev scv 

s.t. ^ + ^ zc/ > 1 v5 c y, 

t;er(S) !7|5c;7 

£ {0, 1} V?; e "l/', 

Z5 G {0, 1} ys c 

We let (LP) be the linear programming relaxation of (|P]), obtained by replacing integrality 
constraints by non-negativity. Its LP dual is as follows. 



max ys + Pv (L)°) 



scv vev 

s ' 



t. Yj ys + Pv<c{v) yvev 

s\veris) 

ucs v£S 

p > 

Call a vertex v £V cheap if c{v) < tt{v), and expensive otherwise. Just like in |T3], we obtain a 



reduced version of (D*^ ) by setting = c{v) whenever w is a cheap vertex, and pv = it{v) otherwise. 
We then define the reduced cost c{v) of vertex t; to be if v is cheap, and we let it be c{v) — Tr{v) 
otherwise. Similarly, we let the reduced penalty Tr{v) be tt{v) — c{v) if v is cheap, and otherwise. 
The reduced dual of Q is then: 

max ^ ys +P{V) (D) 
SCV 

s.t. ys < c{v) yv G V (1) 

s\v&r{S) 

J]yc/<5]vf(^;) ySCV' (2) 

ucs vgS 
y>G 
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1.2 Moss &: Rabani's algorithm 

The algorithm in [14j computes a dual solution for ( |D| ) using a monotone growing process. Initially, 
we let C be the set of all inclusion-wise maximal connected components of the graph induced by 
the cheap vertices and the root. All such components but the one containing the root are active. 

The algorithm raises the dual variables ys corresponding to all active components in C uni- 
formly, maintaining feasibility for ([Pj). When constraint ^ becomes tight for some active com- 
ponent C, C becomes inactive and its dual value is not increased anymore. When constraint ([T]) 
becomes tight for some expensive vertex v, all the (active and inactive) components adjacent to v are 
removed from C and merged to form one new component C given by C := {v}L){C G C : v G T{C)}. 
This new component is active if it does not contain the root, and inactive otherwise. Then, the 
dual growing continues and the algorithm stops when no active component remains. At the end of 
the algorithm, the output is a tree T contained in the connected component in C containing the 
root, and the dual solution {y,p) for (IDj). The claim is that 



with a = O(lgn). 

We exhibit a counterexample where the above process finds a dual solution whose value is a 
factor of ~ n lower than the cost of an optimum solution. The starting point for this example is the 
observation that the set-cover problem is a special case NWST, and that known LP-based O(lnn)- 
analyses for set-cover use dual-fitting, or direct primal rounding approaches. So far, primal-dual 
algorithms based on a monotone dual growing process are only known to have performance ratio 
equal to the maximum frequency of any element (e.g., see |17j). but no better. Since Moss and 
Rabani's algorithm is indeed based on a monotone primal-dual process, this would (interestingly) 
lead to an algorithm of the same type for set-cover. However, as we now show, this is not the case. 

Indeed, the instance shown in Figure [T] draws motivation from the usual set-cover reduction to 
NWST. The graph is obtained by taking a complete bipartite graph with cheap vertices ui, . . . , n„ 
on one side and expensive vertices vi, . . . ,Vn on the other side. For each 1 < i < n, another cheap 
vertex Wi is attached to each expensive vertex Vi. All these cheap vertices have cost and penalty 
n, while all these expensive vertices have penalty and cost n + 1. Finally, we attach to ui an 
expensive vertex of cost 2 and penalty 0, that in turn is attached to the root. 

When running the algorithm of pTlj on this instance, at the beginning the cheap vertices 
ui, . . . ,Un and wi, . . . ,Wn form singleton active components. In the first step, the dual values 
and y{uji} {i = 1, ■ ■ ■ ,n) are increased by 1 and constraint (jT| becomes tight for all expensive vertices 
vi, . . . ,Vn- This basically implies that all such expensive vertices will join together to form a single 
component C containing all vertices Ui,Vi,Wi, for i = 1, . . . , n. The dual value yc of this new active 
component will then be increased by 1 where the expensive vertex adjacent to the root becomes 
tight as well. At this point the entire graph forms one inactive component, and the algorithm 
terminates. The total value of the dual solution y computed by the algorithm is 0{n), while the 
optimal solution is any spanning tree of the whole graph, that has value Q{n'^). 




v<^T 



SCV 
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Figure 1: Counter example for the dual solution proposed in [13j . The pairs of numbers on vertices 
show the values of cost and penalty respectively. 



1.3 The new algorithm 

Our algorithm follows a dual-growing approach, as does that of Moss and Rabani. Unlike their 
algorithm, however, ours is not monotone! Instead, it is inspired by Guha et al.'s primal-dual view 
of Klein and Ravi's NWST algorithm. Like the algorithm in [9J, our method works in phases. In 
each phase but the last we construct trees that merge several of the components from the previous 
phase. When the algorithm terminates, the cost of the tree connected to the root as well as the 
penalty of all vertices not spanned by this tree are approximately charged to a feasible solution to 
(|D^. Both our algorithm and its analysis depart from the standard primal-dual approach taken in 
|9] in several ways. We highlight the two most important new features. 

First, a standard primal-dual strategy for node-weighted problems goes as follows: whenever 
at least two active components tighten constraint ([T]) of a vertex v, then a tree connecting these 
components and v is built. Our algorithm instead considers connecting active and inactive compo- 
nents neighbouring v, but does so only if the total dual of participating components is large enough. 
We stress that the algorithm may build a tree even if only one of the participating components is 
active. 

The second main difference of our methods lies in the charging argument that accounts for the 
cost of the final tree produced. For this, we identify a unique core for each set S in the support 
of the dual solution produced. The cost of the computed final tree is charged to only those cores 
that it spans. We then combine this fact with a non-standard potential function based argument 
to show our claimed approximation guarantee. 

Overall, both our algorithm and its analysis are significantly more involved than the ones of 
Moss and Rabani: in our opinion, this might be consistent with the fact that classical primal-dual 
approaches based on monotone growing seem to fail for NW-PCST, and therefore additional ideas 
are required. 

We describe the details of our algorithm in the next section. We show an execution of our 
algorithm on the example in Figure [T] in Appendix [Aj 
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2 Algorithm 



Our algorithm for NW-PCST is primal-dual, and constructs an integral feasible solution for (|P|) as 



well as a feasible solution for (D^) whose values are within an O(lnn) factor of each other. This is 
accomplished in phases. Phase i starts with a set of initial components as well as a component Tr 
that contains the root. Our algorithm maintains a tree for each of these sets, spanning its vertices. 
We also maintain the invariant, that no two components in are adjacent, or in other words, no 
two components in C* are connected by an edge. 

Initially, in the very first phase of the algorithm we define these sets as follows: recall that we 
call a vertex cheap if its cost is at most its penalty. We look at the graph induced by the root r 
and all the cheap vertices, i.e. G[{r}U{?; : v is cheap}]. Each connected component in this induced 
subgraph will be an initial component in C^, except for the component containing the root that 
will constitute Tr. 



In phase i we run a dual growing process that computes a feasible solution {y\p^) for (D'^), 
and either 

[i] finds a tree T* connecting a set C{T^) C of at least two initial components with cost 
proportional to the total value of the dual solution (y*,^*), or 

[ii] determines that the total penalty of components in C* is at most the value of 

In case [i] we either replace C(T*) in by the single set T*, or we replace by T^- U T*, while in 
case [ii] the algorithm terminates and returns the tree T^. 

To achieve the claimed approximation guarantee, we will show that at least one among the 
feasible dual solutions produced by the algorithm in each phase, has a value within an 

O(lnn) factor the cost of the returned tree Tj. plus the penalties of the nodes \nV \ Tr. 

We now give a detailed description of a phase of our algorithm. For notational convenience 
we will omit superscripts i whenever there is no ambiguity. We will eliminate all p-variables, and 
revert to reduced costs and penalties as described in Section Our dual growing process will 
now compute a dual solution feasible for ([Pj). 

The dual growing procedure of the current phase is best described as a process over time. The 
algorithm maintains a feasible dual solution y'^ for ( [Pj ) for every time r > 0, and we let be a 
collection of sets that contains its support, (i.e., 5"^ contains all S with > 0). 

Call a vertex v tight if the constraint ([T| for v holds with equality for y"^ , and note that cheap 
vertices are always tight. In the following we will call a maximal connected set of tight vertices with 
respect to y'^ a moat. Any two moats are clearly vertex disjoint, and no two vertices in different 
moats can be adjacent. A moat S is active at time r > if 

ucs 

and inactive otherwise. For ease of notation, we let be the collection of active moats. Initially, 
is the set C of initial components, and so each initial component is contained in an active moat 
at time 0. We also define X"^ as the set of all inclusion- wise maximal inactive sets in S'^. For C G C, 
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and T > 0, we let 



age^(C) = min{r, ^ 7r{v)} 



vac 



be the age of C at time r. The age of C is the first time, where C becomes part of an inactive 
moat during the dual growing process if that has happened before time r, and r otherwise. 

For a given set S £ S'^, and an initial component C, we call C the core of /S if C C 5 has 
the largest age value among all initial components contained in S. We will later see that every set 
S £ S'^ has a unique core which we will denote by core(S'). This allows us to extend the age-notion 
to sets S £ S'^: we let the age of S at time r be that of its defining core. 

At any time r > 0, the algorithm grows all sets in uniformly at unit rate. Several events 
may happen during this growth process. 

[A] a constraint of type ^ becomes tight for an active set S S A^, or 

[B] a constraint of type ([T]) becomes tight for a vertex v. 

In [A], the set S now becomes inactive, and moves from to I'^ . If, after S moves, A"^ = 0, 
then the phase ends and we will show that condition [ii] holds. Otherwise, the growing process 
continues for the remaining set of moats in A'^ . 

Consider now an event of type [B]. We say that an initial component C £ C loads v if there is 
a set S" G S'^ with v £ T{S) such that C = core(5). Let C'^{v) be the set of all C S C that load v. 
Uv£ r{Tr) or if 



then we will find a tree T that connects v and the cores that load it, among possibly other things. 
We emphasize the subtle but important departure from the familiar primal-dual theme of active 
mergers: we do not require there to be more than one active moat neighbouring v at time r! 

The tree T is constructed iteratively; initially we let T = {{v},f}>), and we will add to it in 
recursive calls to two main procedures, named FindSubTree (FST) and ConnectVertex (CVtx). We 
will refer to T as the phase tree. 

Once the tree construction is complete, the phase ends. We remove from C all the initial 
components C connected by T and we either replace Tr by Tr U T, in case v £ r(Tr), or add the 
single component T to C. We will show that condition [i] holds for T. 

If instead v ^ r(Tr) and Q does not hold, then we simply continue the growing process. Note 
that, by definition of moat, there will now be a single active moat containing v as well as all the 
sets in A'^ U I'^ adjacent to v. Indeed, at any time r, S'^ is a laminar family of non-adjacent sets. 

We provide more details for the computation of tree T in [B]. The general goal of the compu- 
tation is to construct a tree T with the following property: for every expensive vertex w £ T, T 
connects all the initial components that load wl To this aim, let Si,. . . ,Sp be the collection of 
inclusion-wise maximal sets in that neighbour v. For each of these sets Sj we will now invoke the 
function FST whose job is to find a tree Tj that connects v to those cores in C'^{v) that are contained 
Sj. In general, the procedure FST takes as parameters two vertex sets 5 C y and L C T{S); the set 
L contains vertices that are already spanned by the phase tree T. FST computes a tree connecting 
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L to core(5) and adds it to T. In our computation of tree Tj, we call FST(S'j, {u}). The final tree 
T centered at v will then be the union of {{;} and IJj=i ^i- 

We point out the following subtlety indicated in the figure above. There may be sets S,R £ A'^UZ'^ 
that both load v and hence their cores are contained in C'^{v). Notice however that, if R is contained 
in S, then R is not among the sets Si, ... , Sp, and hence FST will not be called directly for R, but 
rather will be invoked indirectly, at some later point in the recursive procedure. This is important 
to ensure that FST is called at most once for each set in S'^: indeed, this will be crucial for our 
analysis. 

2.1 FindSubTree 

FST(S', L) first constructs auxiliary graph Hs as follows. Let f be the time at which the construction 
of the phase tree started, and hence, where the current phase ended. Let r be the age of S at time 
f . (Note that r could be < r, if S is an inactive set at time r). Start with graph G[S U L] and 
among all the inactive sets in S'^ contained in S, identify inclusion-wise maximal ones; abusing 
notation, we will refer to the super-vertex resulting from identifying such a set -R by i? as well. 

Observe that Hs may now have super-vertices but also original, expensive vertices. In fact, no 
two super-vertices are adjacent, or in other words, every neighbour of a super-vertex in Hs is an 
original expensive vertex. Each such expensive vertex is tight (in terms of constraint ([T])) at time 
f . We let its auxiliary cost cs{v) be the total amount of dual load it feels from core(S'); i.e., 

csiv) = yR- 

RCS,veT(R) 
corc{J?)— corc(5) 

Let the auxiliary costs of the super- vertices and the vertices contained in the core be 0. 
FST has two main parts that we now describe. 

Part I. The purpose of part I is to find a tree Tl that connects the vertices in L, and to add 
it to T. We will later see that L has never more than two vertices. We start by computing a 
minimum-C5-cost path P in Hs connecting the vertices in L. This is trivial when L consists of a 
single vertex a where we let P = {a}. Otherwise, assume that L = {a,b} for some a,b £ r(S'). 
Compute a minimum-c^-cost a,6-path P; i.e., a path in Hs connecting a and b with minimum 
value of 

^csiv). 

We now add the path P to the phase tree T; recall that the two end-points of P are already in 
the phase tree. The path P may contain super- as well as expensive original vertices. Super-vertices 
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(i) (ii) 

Figure 2: Inside method FST. Paths P and Q appear in thick sohd and dashed stroke, respectively. 
Original (expensive) vertices appear in white, and super-vertices are coloured in red; the small red 
subset indicates the core of the considered set. 



will be replaced by sub-trees as follows. For each super-vertex R on this path with neighbours a' 
and b' we recursively call FST with parameters R, and L = {a' , b'}. We also call 

FST(i?',r(i?')nP) 

for each super-vertex R' that is not on P, but is the neighbour of some original expensive vertex 
in P; see Figure [2] for an illustration. The trees returned by all of these subcalls are added to 
the phase tree T, and indeed these trees together with the expensive vertices in P form a tree Tl 
connecting vertices a and b. 

Why do we recursively call the function FST on the sets that neighbour expensive vertices in 
P? The reason is that their cores are loading such expensive vertices. Since those vertices will 
become part of the tree T, as outlined in the previous subsection, we now want to connect all cores 
that load them as well. 

Part II. The goal in the second part of FST is to connect the core of S to the phase tree T. 
This is established by procedure CVtx that takes as parameters a set S, a vertex a G r(S') U S, 
and (for analysis purposes) a level index i. The procedure computes a tree contained in 5 U a that 
connects a to T, and adds that to the phase tree. In our specific case, we choose an arbitrary vertex 
z in core(S'), and call 

CVtx(S', z,0). 

Indeed, at the end, FST adds core(S') to the phase tree as well. 
2.2 Connect Vertex 

Much like FST, CVtx{Sd, Zd,d) first computes an auxiliary graph Hs^, as follows. First, consider 
the graph G[Sd U U T], where T is the current phase tree. Then, as in FST, among all sets in S'^ 
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contained in Sd, identify all inclusion-wise maximal inactive ones, where r = age'^(S'd) and r is the 
time at which the current phase ended. Finally, define the auxiliary cost c^^ (v) of a vertex v to be 
if u is a super- vertex or if v is in core(5d). Otherwise, let cs^{v) be the total amount of dual load 
it feels from core(5rf). 

The procedure computes a shortest-C5^-cost path Qd in Hs^^, from Zd to T. Qd ends in some 
(super- or non-super-) vertex 5^+1 in Hs^^. Note that Sd+i could be a super- vertex. However, this 
would mean that 5^+1 H T 7^ 0, and therefore that the function FST has been called already for the 
set Sd+i- For this reason, we do not call the function FST again on Sd+i, but we rather recursively 
call 

CVtx(5d+i,Zrf+i,d+ 1) 

where Zd+i is the second-last vertex on Qd- Now add the path Qd to the phase tree T. This 
path may contain super- as well as expensive vertices. For all super- vertices S' on Qd we now call 
FST(5',r(5') n Qd)- The computed tree will replace the super- vertex placeholder in phase tree T. 




Each expensive vertex v on Qd may feel dual load from (the inactive component of a) super-vertex 
S", and hence core(5") E J~-^{v)- We observe, however, that among the internal vertices, only Zd+i 
may feel dual load from super- vertices in T n S^. This follows from the fact that Qd is a shortest 
path according to auxiliary costs 05^. The figure above illustrates the situation. CVtx now calls 

FST{S",T{S")nQd) 
for all super- vertices S" not in T that neighbour a vertex v in Qd- 

3 Analysis 
3.1 Correctness 

The goal here is to show that the algorithm described in the previous section can be correctly 
implemented. 

First of all, note that in each phase i the number of initial components in C* decreases by at 
least one, and therefore the algorithm terminates after at most n iterations. 

The final output is a connected tree T^, containing the root by definition. Connectivity of Tr 
easily follows by the way FST and CVtx are defined. However, correctness of the above functions 
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crucially relies on two properties that we just stated in the previous sections and we are now going 
to prove: 

(a) Each set S £ S'^ has a unique core; 

(b) In each call of the function FST the cardinality of the set L is at most 2. 

To see (a), we show the following more general statement by induction: at any time r, an active 
set S G A'^ contains exactly one initial component C with age'^(C) = r. This implies (a), because 
any set in S'^ was an active set up to some time t' < t . Observe that the statement is certainly 
true at the beginning of each phase (r = 0). The only event that changes the family A'^ at time r 
is when a vertex v becomes tight. In order for this to happen, it is necessary that at least one set 
S, with V G r(5), is in A'^ . Suppose we have two such sets Si and 52- By induction hypothesis, 
there is one initial component in Si and one initial component in 5*2 with age r. But in this case, 
condition ([★]) holds and the phase ends. Therefore, if the phase does not terminate, there is only 
one set S G A'^ with v G r(S'). Any other maximal set S' G S'^ adjacent to v, if any, belongs to 
X"^ and by induction hypothesis it follows that it contains only initial components with age < r. 
This implies that the new active set, which now includes S, v, and all its adjacent inactive sets, will 
contain only one initial component with maximum age value r. 

To see (b), we observe that there are only three different situations that lead to FST calls. The 
first one is when Q is verified: here we will call FST with L = {v}. 

The second possibility is when there is a super-vertex R that belongs to either a path P 
computed inside a FST call, or a path Q computed inside a CVtx call. In both cases, we call FST 
with L = {a, b}, where a and b are the expensive vertices on the path adjacent to R. 

The last possibility is when there is a super-vertex R' adjacent to an expensive vertex v that 
belongs to either a path P computed inside a FST call, or a path Q computed inside a CVtx call. 
Assume v £ P (the other case is identical). We call FST with L = {T(R') n P}. What does set 
L look like? Recall that P is a minimum-c^-cost path in an auxiliary graph Hs (for some set 
S). Since super- vertices in Hs are never adjacent, L contains only expensive vertices. We now 
claim that any expensive vertex w G Hs has a positive auxiliary cost cs{w). If not, it means w 
became tight at some time < r' < age'^(S'), where f is the time at which the current phase 
ended. However, when that happened we had to check condition Q: since the phase did not end, 
w became part of an active set at time r', and therefore it would now be inside a super- vertex, a 
contradiction. It follows that, \L\ > 3 then the path P contains at least 3 expensive vertices with 
positive auxiliary cost that are adjacent to R' . Since the auxiliary cost of R' is instead 0, P could 
be shortcut contradicting the fact that it is of minimum cost. 

3.2 Bounding the tree cost in phase i 

Recall that with the exception of the very last phase, each phase i of the algorithm computes a tree 
T* that joins a number of initial components C(T*) in C*. The tree will then be merged with Tr 
if it and have common adjacent vertices, and T* will be added to C*^^ otherwise, replacing the 
components of C(T*). In this section we will bound the total cost of T*. We do this by providing a 
detailed charging scheme that distributes c(T*) over the components in C(T*). 
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Recall from the description of phase i in Section [2] that we start construction of tree T* as soon 
as constraint ([T]) becomes tight for some vertex v and either v G ^{Tr) or ([★]) holds. Suppose that 
this happens at time r > 0, and recall that y'^ is the feasible dual solution for (|D]) at this time. 
By definition the set contains all sets S in the support of y'^, and it is not hard to see that the 
family of these sets is laminar (i.e., any two sets are either disjoint, or one is fully contained in the 
other) . 

As soon as a vertex v becomes tight and either v G r(Tr) or Q holds, our algorithm invokes 
FST for each inclusion- wise maximal set S € S'^ that neighbours v; the goal being to connect the 
core of S to v. Each of these top-level FST calls may itself trigger further, lower level FST calls on 
sets S' that are descendants of S in the laminar family defined by 5"^. 

In Part I, Function FST constructs a partial tree to be added to T. In part II, FST invokes 
CVtx in order to connect vertex zq G core(S') to the partial tree Tl constructed so far. In its 
construction of a path linking core and T^, CVtx will call itself with progressively higher level- 
indices. A chain of CVtx calls is associated with a chain {/SjI^^q of sets in S"^ with the following 
properties: 

(i) So = S, 

(ii) Si is a super- vertex in Hsi_^ contained in Si-i for all 1 < i < p — \, 

(iii) Sp is an original vertex contained in Hsp_^- 




Figure 3: The chain of CVtx calls invoked from within a call to FST. The thick black lines indicate 
edges of the tree Tl constucted by FST. 
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For < i < p — 2, the level-i cah to CVtx computes a shortest C5.-cost Zi, Sj+i-path Qi in Hs^. 
By assumption Si+i is a super- vertex of Hs^- Let Zj+i be its original, expensive predecessor on 
path Qi, and let Qi denote its Zj, Zj+i-prefix. We then inductively call 

CVtx(S'j+i,Zi+i,i + 1). 

Finally, Sp is an original vertex that is part of tree Tl, and the zq, T-path is simply the concatenation 
of Qoi • • • J Qp-i- Figure [3] illustrates the construction for p = 3. 

A component C G C* is charged whenever FST is called on set S whose core is C. As we will 
see later, during a phase, FST is executed at most once for any given set S, and thus, a component 
C £ C is charged at most once during a phase i. Even stronger, we will see that any core C is 
charged at most once throughout the entire execution of the algorithm. This fact will be crucial 
for the performance ratio analysis in the next section. 

As mentioned, a call of the form FST(5, L) yields a charge 'I>(core(S')) to core(5'). This charge 
consists of three main components. 

(CI) The auxiliary cost cs{P) of the L-path P computed. 
(C2) The auxiliary cost 

j=0 

of the paths computed in the CVtx calls. 
(C3) The non-auxiliary costs 

p-i 

^ic{zj) - Cs^_^{Zj)) 
3=1 

of vertices zi, . . . , Zp-i. 

Several comments are in order. Consider an expensive vertex t; on P U Qo U . . . Qp-i. We charge 
the vertex' auxiliary cost to core(S'). Who pays for the non-auxiliary part of v's cost? Recall that, 
in the description of FST and CVtx, v 7^ zi, . . . , Zp_i we are careful to invoke FST for all super- 
vertices S' whose cores are in C^{v). The set C^{v) contains all cores of sets in 5'^ that ever load v. 
Thus, through these FST calls, we make sure that the entire cost of all vertices in the constructed 
trees is charged to some cores. 

Vertices zi, . . . , Zp-i receive special treatment; why? A vertex Zj, 1 < j < p — 1 \s incident to 
super-vertex Sj in Hsj-i where Sj intersects T. Vertex zj may also be adjacent to other super- 
vertices that intersect T. All of these super-vertices have cg^.j-cost 0. If a vertex Zj feels dual load 
from super vertex R that intersects T, then we may have already called FST on R. We must not call 
FST twice for a set i? G 5"^! For this reason, we charge all its reduced cost c{v) to core(S') instead. 
Differently, note that Sp is an expensive vertex by definition, and the vertex Zp preceeding Sp on 
Qp-i is not incident to super-vertices intersecting T. Otherwise, Qp-i would not be a shortest 
Zp-i, T-path in i?5p-i- 

We now bound the total charge $(C) to a core C in phase i. Let us first show that each core 
is charged at most once in the entire algorithm. 
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Lemma 3.1. Any set C of vertices is charged at most once in the entire algorithm; i.e., FST is 
called at most once on a set S whose core is C . 



Proof. First note that whenever FST is called for a set S with core C in some phase, then C is 
connected to at least one other core through the phase tree. Hence, C will be removed from the 
list of initial components C in this phase, and will never enter it again. Hence, it suffices to show 
that a set C G C is charged at most once in any given phase i. 

Within phase i, let S be the first set with core C for which FST is called. FST(5', L) recursively 
computes a tree within G[SyjL\ that connects the core of S (i.e. C) to the phase tree. All FST calls 
from within FST(S', L) must be on inactive strict subsets of S, and hence all these subsets have a 
different core. On the other hand, after FST(5', L) is completed, core(S) is connected to the phase 
tree. 

Finally, note that in our algorithm, we never call FST for sets S whose core is part of the phase 
tree. This completes the argument. □ 

We continue by focusing on a specific component C G C*. Suppose that C was charged during 
FST(S', L). Recall that we let r be the time at which the construction of the phase tree started, and 
hence, where the phase ended. We will bound the part of $(C) coming from (CI), (C2), and (C3). 

Lemma 3.2. The auxiliary cost cs{P) of the path connecting the vertices in L is at most 2 age^{S). 

Proof. By the standard primal-dual argument on shortest paths, it follows that any path from any 
vertex in Hs to any vertex in core(S') has auxiliary cost at most age'^(S'). Therefore, the shortest 
path between any 2 vertices in Hg has auxiliary cost at most 2age'^(5'). □ 

In order to bound (C2), we use a similar primal-dual argument, as well as Q. 

Lemma 3.3. Suppose that the path constructed in the CVtx call to connect C to the phase tree 
consists of the concatenation of Qo, . . . , Qp-i. The auxiliary cost Y2^=o (^Sj (Qj) of this path is at 
most 4:age^{S). 

Proof. Let 5o = S, Si, . . . , Sp be the chain of sets in 5"^ corresponding to the chain of CVtx calls 
invoked to construct Qq, . . . , Qp-i. Using a similar argument as in Lemma 3.2 one sees that that 



C5,(Q.) <2age^(5i)- 

The lemma follows by showing that age'^(S'j) < age'^(S'i_i)/2 for all 1 < i < p — 1. 

Observe that Si is a super-vertex in Hsi_^ for all 1 < i < p — 1. Hence, there is a time r < r 
during the algorithm where the active moat containing C meets Si. Repeating an argument used 
in Section 3.1 , Si could not have been active at the time as Q would have been satisfied otherwise. 



and the phase would have ended. More specifically, the age of Si must have been smaller than 
f/2 < age"(5,-i)/2. □ 

It remains to bound the contributions in (C3). 
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Lemma 3.4. Using the notation of Lemma 3.3, let Zj+i be the final vertex of path Qi for all 
< i < p — 1. The total non-auxiliary cost 



p-i 

i=l 

is at most age'^[S). 

Proof. This once again uses condition Q. Since is a (super-)vertex in Hs^ for all < i < p — l, 
it must be the case that at some time f < age'^(S'j), an active moat containing core(5j) tightens 
vertex Zj+i (i.e., constraint ([T]) of Zj+i becomes tight). Once again, as the phase did not end 
at this time, we know that the total age of components with cores in /^"^(zj+i) — core(S'i) (and 
hence the non-auxiliary reduced cost of the vertex) must be less than f/2 which in turn is at most 
age'^(S'j)/2. As age'^(S'i) is at most age'^(S'o)/2*, and hence the non-auxiliary costs of vertices Zi 
decrease geometrically, the lemma follows. □ 

We obtain the following final corollary. 

Corollary 3.5. The total charge $(C) of a component C £ is at most 7 age'^{C). 

Note that, trivially, age'^(C) = YliS-c=coie{S) Vs- Finally note that the cost of each vertex in 
the phase tree is charged to some core in C*. Hence, the reduced cost of the T* is at most the 
sum of <I?(C) over all initial components C C C{T^). 

3.3 Approximation Factor Guarantee 

In this section we prove that our algorithm is a primal-dual O(lnn) approximation algorithm for 
the Prize Collecting Steiner Tree problem. 

We start by proving that the total reduced cost of the final tree T returned by the algorithm in 
the last phase is within an O(lnn) factor of the value of a feasible dual solution y* to ( [Pj ) computed 
during some phase of the algorithm. 

Indeed, for the purpose of designing an LMP algorithm the following section we will prove a 
slightly stronger statement. The algorithm computes a feasible dual solution for every phase 
i = 1, . . . , m where m is the total number of phases. For any such t/*, we consider the dual solution 
y*, obtained from y by setting y^ := ys if core(5) C T, and yg := otherwise. Clearly, is a 
feasible dual solution. 

We will prove the claimed bound by showing that it is possible to select / = O(lnn) indices 
ti, . . . , E {1, . . . , m} such that 

-<T) = o(j2Y.ys)- (3) 

i=i s 

From the previous section, we know that 

-c{T) = -c{r) < ^(^) 

i-.T^QT C&C^:CCT 
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Therefore, it will be enough for us to prove that 

C£Ci:CCT j=l S 

To this aim, let us call C* is the collection of initial components in phase i that have been 
included in the final tree T, that is C* = {C £ : C CI T}. 

The proof will go as follows: first, we will define buckets 1, . . . ,1 and assign every component 
in Ui^* to exactly one bucket. Secondly, we will prove that for every bucket j, we can identify 
one dual solution y*^ such that the total charge of the components in bucket j is within a constant 
factor of the value of y*-' . 

Let us start by describing the assignment process of components to buckets. 

Set I := [InnJ. In order to avoid dealing with constant in the formulas later on, let us set 
$'(C) := ^{C)/7 for every C e Ui^* and let ^'^ax be the maximum <i>'(C) among ah C G [JiC\ 
We assign a component C to a bucket 1 < j < / if and only if 

max , r^s'/'/^^ ^ max 

and we assign C to bucket / otherwise. Clearly every component is assigned to exactly one 
bucket. We observe that bucket 1 is non-empty by definition, and bucket I contains all small 
charges. 

We now describe how we select one dual solution for each bucket 1 < j < /. Let tj be equal to 
the smallest phase index i such that there is some component C G C* assigned to bucket j. If there 
are buckets with no components assigned to it, then we let tj be any arbitrary phase index. Note 
that for some j and j' ^ j, we may have tj = tji. We also let ti = ti, and we will show that 

y: i>'(c)=o(5]y|), (o) 

C assigned to ^ 
bucket j 

for aU l<j<l. 

From the discussion above, it follows that Q implies the bound ([s]) on c(T). Therefore our 

goal now is to prove the above equality. We will treat the cases j < I and j = I separately. 

<j)' 

Let us start by assuming j < I, and let bj := — Observe that, by construction, every 
component C assigned to bucket j satisfies: ^>'(C) < 2bj. Therefore, if bucket j contains in total 
K components, the left-hand side of Q is at most 2Kbj. The key idea to relate the quantity 2Kbj 
to the value of the dual solution y*J is that of introducing a potential function /3j{i), defined as 

/3,(i)= ^ min{7f(C),6,}, 

for every bucket j, and for each phase i. 

We now show how to use this potential function. Before, we just state a useful remark that 
follows directly from Corollary |3.5| and the definitions. 
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Remark 3.6. Let r* be the time when phase i terminates. If there is a component C G C which 
is assigned to bucket j , then f* > age^\C) > <I>'(C) > hj. Moreover, 7f(C) > <I>'(C). 

The next lemma shows that the value of the potential function computed in phase tj is a lower 
bound on the value of the dual solution y*J . 

Lemma 3.7. f3j{tj) < YlsV^s- 



Proof. By Remark 3.6 above, the time r := r*-' when phase tj terminates is at least bj. So for every 
component C £ C^^ we have that the age of C at time r is equal to min{7f(C), r} > min{7f(C), 6j} 
and therefore, 

Y.ys> E age^(C)> E min{6„7f(C)} = /3,(t,). 



□ 



The following lemma is the heart of our analysis. It shows that in each phase i, the potential 
function /3j decreases by an amount proportional to the total charge value of the components in C* 
assigned to bucket j. 



Lemma 3.8. Consider a phase i < m, and let k be the number of components in C assigned to 
bucket j . Then 

k, 
2 



/3,-(i)-/3,(2 + l)> 



Proof. Let r := r* and T* be the subtree constructed in phase i. Moreover, let C{T^) = {C G C : 
C C T*} be the set of initial components of phase i connected by T*. If C{T^) = 0, then C* = C*^^ 
and /c = 0, therefore the statement holds. 

So, suppose C(T*) ^ 0. Note that the only difference between C and C is that C does 
not contain the components in C{T^). Instead, either C contains the new component T* or T* 
joins Tj.. Therefore in both cases we have 

/5,(i)-/3,(i + l) = 



^ min{6j,7r(C)} - ^ min{6j, 7r(C)} 



> mm{bj,7t{C)} -mm{bj,Tr{T)}. 

C£C{Ti) 

Let the value of the right hand side of the equality above be X. We distinguish 3 cases depending 
on the value k of components in assigned to bucket j. 

If /c = 0, it suffices to show that Pj{i) is non-increasing in i. This is immediate as 

X > mm{bj, 7f (C)} - min{6j, vf (T*)} = 0. 

C6C(T») 
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If A; > 1, let C G be any component assigned to bucket j. Then, it{T'^) > 7r(C) > <I>'(C) > bj by 
Remark 3.6 and therefore: 

x>{ 

C£C{Tiy.^'(C)>bj 

k 

> kbj - bj >{k- l)bj > i^bj. 

If A; = 1, again vf(T*) > bj. Since k = 1, there is exactly one component (7 G C* merging into 
T* that is assigned to bucket j, and we have vf(C') > ^'{C) > bj. Moreover, for every C G C{T^) 
we have m.m{Tt{C),bj} > min{age'^(C), 6^} and since the termination condition Q is satisfied, we 
have: 

X > (min{7f(C), bj} + ^ m.m{Tt{C),bj}) - bj 

C&C{T^)-C 

> {bj + min{age"(C), - bj 

> min{ Y age^(C),6j} 

= min{ J2 age^(C) - age^(C'), 6^} 
CeC{T') 

3 Ik 

> min{- r - T, bj} > - bj > - bj 

□ 



With the above two lemmas at hand, we are now ready to prove Q. Let ki be the number of 
the components in C* assigned to bucket j. For a bucket j we have: 



m—l . 

Yy'^>m)>Y^bj + ^,{m) 

S i=tj 

C assigned to 
bucket j 

and hence 

Y ^'{c) = o{Yy's)- 

C assigned to 5cy 
bucket j 

It remains to prove Q for j = I. At the be ginning of the algorithm, we have at most n initial 
components, and in every phase, at least two of these are merged into one common component. 
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Hence, the total number of distinct initial components throughout the algorithm is at most 2n, and 
this is an upper bound on the total number of components that are charged in the algorithm. 

We know that in phase ti = ti there is a component C G C '' such that <I>'(C") > <l>^^2,/2. Let 
T := f*^ Then, Es^i > age-(C') > <!>'{C') = a>'^„j2. Therefore 

C:'S>'{C)<bi S 

as desired. Putting all together, we proved: 

Theorem 3.9. Let T be the tree returned by the algorithm and let y* be the dual solution to (^ | 
among y^, ■ ■ ■ ,|7™ with maximum value. Then 



c{T)<0{lnn)Y,y*s- 



The next lemma will complete our argument. Its proof follows trivially by the algorithm 
definition. 

Lemma 3.10. Let y"^ be the dual solution found in the last phase of the algorithm, then every 
component S in the support of this solution is disjoint from T and 

vf(F'\T)<^y^ 
s 

The approximation bound now follows. 

Theorem 3.11. Let OPT be the value (cost plus penalty) of the optimal solution for Prize Col- 
lecting Steiner Tree problem on an instance of problem with n vertices, then, Algorithm^^ finds a 
solution T with c{T) + Tr{V' \T) = 0{liin)OPT. 

Proof Recall that, given a feasible solution y to the dual ([D|, setting p^ = c{v) if v is cheap, and 



Py = 7r{v) otherwise, yields a feasible solution {y,p) to the dual (D''). Using weak duality together 
with Theorem 3.9 and Lemma |3.10[ i we have 



c(r) + Tr{V' \T)<Y,Pv + c{T) + 7f,(y' \ T) 
= 0(ln n)OPT + OPT = 0(ln n)OPT 

□ 

4 LMP Algorithm 

A Lagrangean Multiplier Preserving (LMP) a-approximation algorithm is an algorithm that finds 
a solution F for an instance of the problem such that 

c(F) + aTT{V \F)< aOPT, 
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where OPT is the value of the optimal solution for the considered instance. In the following, we 
show how the algorithm described in the Section [2] can be used in a black-box fashion to obtain an 
LMP 0(lnn)-approximation algorithm for NW-PCST. 

Theorem 4.1. There is an LMP O (Inn) -approximation algorithm for NW-PCST. 

Proof. Assuming that Algorithm [T] is an a- approximation algorithm for NW-PCST problem, we 
will show how to obtain an LMP 2Q!-approximation. 

Given an instance {G, tt, c) of NW-PCST problem, consider the following algorithm: 

1 Set tt'{v) = 2Tr{v) — c{v) for every cheap vertex v and it'{v) = 7r{v) for every expensive vertex 
v; 

2 Run Algorithm [l] for the instance (G, 7r',c); 

3 Output the tree T returned by Algorithm [T] . 

Observe that if we denote the reduced penalties of the new instance by vf', then Tt'{v) = 2Tt{v) 
for every vertex. 

By Theorem 3.9, we can find a dual solution y to ( [Pj ) for the instance (G, vr', c), constructed in 
some phase i < m of the algorithm, such that 

c(T)<a (4) 

SCV':cove(S)CT 

Furthermore, if z is the dual solution found in the last phase of the algorithm, then by Lemma 
13.101 we have 

7r'{V'\T) = 27r{V'\T)< J] zs- (5) 

scv 



Claim. The vector y' obtained by setting y's '■= ^ + ^ for all S C V' , is a feasible solution 
to ( [D| ) for the original instance (G, 7r,c). 

It is easy to see that the cost constraints ([T]) are satisfied, since the costs (and therefore the 
reduced costs) never change and y' is simply a convex combination of two solutions satisfying that 
set of constraints. 

So, we only need to prove that the penalty constraints ^ are satisfied for every set S Q V' . 
Let C be the set of the initial components considered by the algorithm at the beginning of the phase 
i (where solution y was constructed), and let T> be the set of initial components considered by the 
algorithm at the beginning of the last phase (where z was constructed). Recall that each set in 
the support of y (resp. z), has a single core, and the age of the set is not greater than the total 
reduced penalty of its core. Moreover, by construction, the core of every set in the support of y is 
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contained in T. Then, it follows that 



^ + zr/2) = Vr/^ + E 

RCS RCS RCS 



C&C:CCSnT DeV:DCS 

= E *(^)+ E *(^) 

ceCiCcsnT Dev-.Dcs 

Observe that z is the dual solution constructed by algorithm during the last phase m, and therefore 
any initial component in T> is either an initial component in C as well, or it is the union of a bunch 
of initial components in C. Since all components in T) became inactive, none of them intersects T. 
Therefore Edsb^dcs ^(^) = Ec6C:Ccs\T ^(^)- We obtain 

Ey«< E ^{C) = ^{S). 

RCS C&C:CCS 

Therefore, the penalty constraint holds for S and the claim is proved. 
Adding a times the inequality ^ to the inequality Q, we get 

c{T)+2a7t{V' \T) <aJ2ys + aJ2 

scv scv 

scv scv scv 

and hence 

c(r) + 2an{V' \T) < 

scv vev 

which means the algorithm is an LMP 2a-approximation algorithm. The last inequality follows 
again from the fact that setting = c{v) if v is cheap, and py = 7r{v) otherwise, yields a feasible 
solution to the dual (D°). □ 
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A Example of Execution 



We apply our algorithm to the instance shown in Figure [T} 

Phase 1. We have = {r}, while each cheap vertex Uj and Wj, j = l,...,n, form an 
active (singleton) component {uj} and {wj}, that is in . We raise the dual variables of the 
active components till, at time r = 1, we check condition Q for the expensive vertex vi that got 
tight (indeed, since all the expensive vertices Vj are tight, we could arbitrarily pick one of them). 
Condition Q holds, and the tree we construct is simply the star centered at vi. At this point 
the phase ends, and note that our first dual solution is exactly the dual solution found by the 
algorithm of fL3] . 

Phase 2. We still have Tj. = {r}, while is now different. contains the active component 
spanned by T^, as well as an active (singleton) component {wj^ for all j 7^ 1. We raise the dual 
variables of the active components till, at time r = 2, the expensive vertex x adjacent to the root 
becomes tight. The tree we construct is simply U {x}. At this point the phase ends, and 
note that the value of the dual solution is again 0{n). 

Phase 3. We now have Tr = {r} U T'^, while contains an active (singleton) component {wj^ 
for all j ^ 1. We raise the dual variables of the active components till, at time t = n, they all 
become inactive! At this point the algorithm ends, by returning the current Tr- Note that the 
value of the solution found is c{vi) + c(x) + X^j=2 '^{wj) = (n + 1) + 2 + n • (n — 1), and that the 
value of the dual solution is indeed n • (n — 1)! 

B Algorithm Pseudocode 
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Algorithm 1 PrizeCollectingStemer(G(y, £■), c, 7r,r) 

^ {S : 5 is a connected component not containing r 

in the graph indTiccd by cheap vertices and r} 

Tr the connected component containing r in the graph induced by cheap vertices and r. 
Let i denote the phase number, 
while ^ do 

i ^ i + 1 

Initialize y'g ^ for aU S C V', A' ^ C\ T ^ 0, T' ^ , ^ 
while A'^^do 

ei ^ mm^gy\u^,g^.^^.5i |{se^':^er(5)}| > 
{; ■<— the vertex minimizing the statement above 
e2 ^ minsg^,{^^g5 7f(v) - Eiic5 2/il} 
e ^ min{ei, 62} 
r' + e 

2/5 ^ + e for all 5 G 

agc^(S') age^(core(S')) ^ for all MS G and r > 
if e = €2 then 

for5G^^:7f(5) = Eficsy}ido 

Remove S from and add it to I* 

end for 
else 

H ^ {inclusion-wise maximal 5 G .4' U I* : v G r(5)} 
if Ecez:ri(c) age^'(C) < | and {5 ^ r(r,) then 

Remove all sets in from A!' and X' 

add {i;} U ( \Jr^ R) to A' 
else 

for 5 G do 
FST(5,{i}}) 
end for 
break while 
end if 
end if 
end while 
if 7^ then 

^i+i ^c\{C eC -.C CT'} 
if n r(r) / then 

Add vertices of T' to Tr 
else 

Add to C^'^^ as a new initial component 

end if 
else 

return 
end if 
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end while 
return Tr 



Algorithm 2 FST(S', L) 

Construct the auxihary graph Hs and compute the auxihary costs c'g 
if \L\ = 1 then 

else 

Assume that L = {a, b} for some a,b G r(S') 
P a minimum c'g cost a,6-shortest path in Hs 
end if 

for each super vertex R G Hs : P fl r(i?) / do 

FST(i?,Pnr(i?)) 
end for 

■<— U { original vertices in P} 
zq ^ arbitrary vertex in core(iS') 
CVtx(5,zo,0) 
Ti-^TiU core(S') 



Algorithm 3 CVtx(S', z, d) 

Construct the auxiliary graph Hs and auxiliary costs c'g 

For every super vertex R G V{Hs) if fl Tj 7^ 0, identify the vertices of i? fl T to super vertex 
Rinf 

Q = qi . . . qi a, minimum c'g cost ZjT-shortest path in Hs 
Tj Tj U { original vertices in Q} 

for super vertex R G V{Hs) \ V{f) : Q n r(i?) 7^ do 

FST(i?, r{R) n Q) 
end for 

if = i? is a super vertex then 

CVtx(i?,gi_i,d+l) 
end if 
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